Magnetic-field generation in helical turbulence 
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We investigate analytically the amplification of a weak magnetic field in a homogeneous and 
isotropic turbulent flow lacking reflectional symmetry (helical turbulence). We propose that the 
spectral distributions of magnetic energy and magnetic helicity can be found as eigenmodes of a 
self-adjoint, Schrodinger-type system of evolution equations. We argue that large-scale and small- 
scale magnetic fluctuations cannot be effectively separated, and that the conventional a-model is, 
in general, not an adequate description of the large-scale dynamo mechanism. As a consequence, 
the correct numerical modeling of such processes should resolve magnetic fluctuations down to the 
very small, resistive scales. 
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1. Introduction. It is well established both analyti- 
cally and numerically that a weak magnetic field can be 
amplified by the random motions of a highly conducting 
fluid HSIa. This occurs because magnetic field lines are 
generically stretched by the random motions of the fluid 
in which they are (almost) "frozen" . Such mechanisms 
of turbulent dynamo action are invoked to explain the 
origin of magnetic fields in astrophysical systems, such 
as planets, stars, the interstellar and the intergalactic 
medium, etc. 

In many cases magnetic fields are observed to be strong 
and ordered on scales much larger than the velocity cor- 
relation length. The traditional view being that the ori- 
gin of these fields can still be explained in the frame- 
work of isotropic and homogeneous turbulence provided 
the latter lacks reflectional symmetry. For this case the 
helicity integral of the velocity can be non-zero, i.e., 
H = J v • (V x v)dV ^ 0. 

To illustrate this idea let us assume that the magnetic 
field, compared to the velocity correlation length, l ve i, 
possesses only large scale and small scale components. 
The magnetic field evolution is described by the induction 
equation 



d t B = V x (v x B) + 77 AB, 



(1) 



where rj is the (collisional) diffusivity. Averaging over the 
small-scale fluctuations, I < l ve i we obtain the equation 
for the large-scale, or the mean magnetic field B(a;, t). It 
can be written in the general form 



9 f B = Vx (aB) + /JAB, 



(2) 



where it is assumed that the mean field varies slowly in 
space, and, therefore, its higher-order spatial derivatives 
can be neglected. The parameters a and (3 can be esti- 
mated on dimensional grounds to be a ~ (v • (V x v))r„, 
and j3 ~ (v 2 )r„, where t v is the velocity correlation time. 

In Fourier space, the linear equation @ has the eigen- 
values, Ai = —fik 2 , and A2 3 = —/3k 2 ± ak, where k is 



the wavenumber. Thus a growing eigenmode always ex- 
ists provided small enough wave numbers are allowed in 
the system under consideration (galaxy, laboratory de- 
vice, simulation box). The maximal growth rate is then 
given by 70 = a 2 / (4/3), and the corresponding scale of 
the growing mean magnetic field is Iq ~ 2(3 / a. In order 
to comply with the underlying assumption of scale sep- 
aration it is assumed that this scale is much larger than 
the velocity correlation scale l ve i. 

The mean-field growth rate, 70 vanishes if the velocity 
fluctuations possess no helicity, 7i = 0. One therefore 
might expect that the generation of magnetic fields at 
large scales, I > l ve i, in homogeneous and isotropic tur- 
bulence may only be possible if the velocity field lacks 
reflectional symmetry, and that such magnetic fields are 
described by the mean-field equation J2J). This effect is 
traditionally called the a-dynamo mechanism. 

However, numerical results suggest that the large- 
scale magnetic field evolution in helical turbulence may 
not be adequately described by the a- mechanism iffi). 
For example, Vainshtein and Cattaneo Q noted that 
the small-scale magnetic fields are amplified more effec- 
tively than the large-scale ones, and when their energy 
is large enough to affect the velocity dynamics, the a- 
mechanism may become much less effective. The influ- 
ence of small-scale magnetic fields on the large-scale dy- 
namo mechanism @ has been stressed in many works 
(see, e.g., jHEl). 

Previous investigations of inconsistencies related to the 
a-dynamo mechanism essentially concentrated on the 
nonlinear effects related to dynamo saturation, which, so 
far, have resisted exact analytical treatment. Present- 
day direct numerical simulations of turbulent dynamo 
action cannot provide conclusive results either, due to 
quite limited numerical resolution, e.g., [llj . 

In the present paper we propose that some of the es- 
sential physics of large-scale magnetic field generation is, 
in fact, captured already at the initial, kinematic stage of 
dynamo action. Our analysis is based on the exactly solv- 
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able model of dynamo action due to Kazantsev • De- 
spite the many simplifying assumptions about the statis- 
tics of the velocity field this model has proven to be a 
valuable tool in understanding the dynamo mechanism. 
In particular, it treats the induction equation exactly, 
and it allows a rigorous derivation of the a-model equa- 
tion (J2J) . So far only the non-helical case has been exten- 
sively analyzed in the literature, e.g., 0,0,0]- Here 
we address the problem in its full generality. 

As an important new result, we show that the evolu- 
tion equations for the magnetic-energy and the magnetic- 
helicity have self-adjoint structure (we note here that 
although the kinematic-dynamo equations have been 
known for over 20 years, their self-adjoint structure had 
not so far been discovered). As a consequence, in the 
kinematic regime, the spectrum of magnetic fluctua- 
tions can be expressed as a sum of eigenf unctions of a 
Schrodinger-type equation with imaginary time, where 
the eigenvalue, A, gives the growth rate of the correspond- 
ing mode. 

In analogy to the quantum-mechanical states in a po- 
tential well, the eigenmodes growing with A < 270 corre- 
spond to "travelling particles," i.e. they are correlated at 
the system size. By contrast, the faster-growing modes 
(with A > 270) correspond to "trapped particles" ; their 
correlation lengths are less than infinity, and they fill the 
whole range of scales down to the resistive scale. At any 
given scale the modes with A > 270 may rapidly become 
dominant over the slowly growing non-localized modes. 

The eigenmodes with A > 270 are not captured by 
the mean-field equation (J2J), consequently the a-dynamo 
model P)l. based on the assumption of scale separation 
and on small-scale smoothing is inadequate. To describe 
the large-scale dynamo mechanism correctly, numerical 
simulations of uniform, isotropic, helical turbulence must 
resolve the full range of scales from to to l v . Further- 
more, the origin of large-scale fields, such as those ob- 
served in astrophysical situations, may be related to the 
nonzero large-scale average of the fluctuating part of the 
field, and not to the mean-field as described by mean-field 
models. 

2. The Kazantsev model for helical kinematic dynamo. 
Kazantsev and Kraichnan 0] introduced the solv- 
able models in the theory of passive random advection. 
The essential assumption is that the random velocity field 
is Gaussian and short-time correlated. It is also assumed 
that the velocity field has zero mean, (v) = 0, so that 
the problem is completely specified by the velocity co- 
variance tensor. For the statistically homogeneous and 
isotropic case, the covariance can be written as 



(t/(x,iy(x',f')> = K ij (\x-x'\)5(t-t'), 



(3) 



where k %3 is an isotropic tensor. For mirror-symmetric 
velocities, the correlation tensor n lJ is symmetric with 
respect to the interchange of the indices i and j. In the 



general case, however, this tensor has both symmetric 
and antisymmetric parts, 



k 13 {x) =k n [5 13 



KL- 



ge J x 



The first two terms at the right-hand side of I0J repre- 
sent the symmetric, non-helical part, while the function 
g (x) describes the helical part of the velocity correla- 
tion tensor. Here e l3k is the completely anti-symmetric 
pseudo-tensor, and summation over the repeated indices 
is assumed. The requirement that the velocity be incom- 
pressible implies that kn(x) = kl(x) + xk' l (x)/2, where 
the primes denotes derivatives with respect to x. 

The magnetic field correlator can similarly be intro- 
duced: H i3 (x,t) = (5*(x,i).BJ(0,t)), satisfying 



H 13 = M N ( S 13 



M L - 



XX- 



Ke ijk x k , 



(5) 



where the corresponding solcnoidality constraint implies 
A/jv = Ml + xM' L /2. Our goal is to find the functions 
Ml(x, t) and K(x, t) that contain the information about 
the magnetic energy and the magnetic helicity. 

Differentiating H 13 (x,t) with respect to t and making 
use of (JTJ, ©, and we obtain after a cumbersome 
but straightforward calculation, that the magnetic corre- 
lation tensor obeys 



d t H i3 = Jl imn Jl3 rt ^pmrjjnt\ 



(6) 



where R v 



= e lW e imn Vfc, and = d/dx k . An analo- 
gous, although not identical, representation of this equa- 
tion was derived in [J^, but see also the derivation 
in 0, 0] . The tensor T 1J can be represented in the fol- 
lowing form: 



rplj 



where 



A ( 5 i3 

V2 V 



A(x) 
B(x) 
C(x) 



x l x 3 



+ B 



XX 



C 



— Kif(x) + 2r)) , 
{kl(0)-k l (x) + 2t)), 
V2(g(0)-g(x))x, 



(J) 



(8) 
(9) 

(10) 



and symbol B(x) in © should not be confused with the 
magnetic field B l (x,t) in Eq JSJ. Hereinafter, we adopt 
the notation kq = kl(0) = kjv(0), and go = g(0). 

Equation © can be considerably simplified, since the 
magnetic-field tensor (JjJJ contains only two independent 
functions, Mr, and K. The reduced equations were de- 
rived in |l4|. however, the symmetric structure of the 
tensor equation © was not preserved. In the next sec- 
tion, we derive the reduced equations, keeping their sym- 
metric structure intact. In this way, we reveal the self- 
adjoint nature of the equations which allows us to gain 
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new insight into the large-scale dynamo mechanism, and 
to elucidate the limitations of the conventional a-dynamo 
paradigm presented in the introduction. 

3. The self-adjoint dynamo equations. In this section 
we show that (0 is self-adjoint. We begin by rewriting 
(|rj|l in the equivalent form, 



dt&i = DpgH nt , 
where D is the self-adjoint differential operator 

and the matrix J is symmetric, 

jll npl tqlrppq 

J nt — e e 1 1 



(11) 



(12) 



(13) 



i.e., this matrix does not change under the interchange 
of its lower and upper sets of indices. We now express 
the operators D and J in the basis defined by the three 
orthogonal "vectors" : 



^2 



1 



cV2 



1 e ijs xS 

xV2 x 



(14) 

(15) 
(16) 



This is possible since the functions k u , and can 
themselves be expanded in this basis. The normaliza- 
tion of the vectors £,f = £f = £f = 1/x 2 is chosen in 
such a way as to preserve the self-adjoint structure of 
the differential operator, D, as we will see presently. A 
straightforward calculation leads to: 



J = 



B A 
A C 
C B 



(17) 



\xM' L , and 



{y/2xM N ,xM L ,y/2x 2 K}, where M N = M L 
Ml and K are some independent functions (compare this 
with 10). We further require that the vector W satisfies 
the equation: 



d t W = —R T JRW, 



(20) 



then, clearly, the function H obeys the dynamo equa- 
tions (|lip. as can be verified by applying the operator R 
to both sides of Eq. I|2l)|l . The operator in the right-hand- 
side of l|20(l is now explicitly self-adjoint. This represen- 
tation constitutes the formal solution of our problem. 

For practical purposes, 1)20(1 can be further simplified 
since only two components of the vector H are indepen- 
dent. Conveniently, the necessary reduction is already 
present in i|20|) . Indeed, calculating the operator in the 
right-hand-side of l(20|l . one sees that it acts only on the 
second and the third components of the vector W, so the 
system is automatically reduced to the two independent 
equations that preserve the initial symmetry structure. 
The reduced equations have the self-adjoint form 



d t W= —R JEW, 



(21) 



where E is the reduced form of the operator E, and J is 
the reduced form of J: 



R = 



X 









j = 



E C 
C B 



Here we introduced the self-adjoint operator 

xA'). 



E = 



1 9 o 9 1 /, 

2 dx dx y/2 



(22) 



(23) 



The validity of (|21l I22[) can be verified most easily by 
direct calculation of the right-hand sides of (|2U|I and 12111 . 
For convenience, we write out the matrix form of Eq. 121|) 
explicitly: 



D 



d 2 

V2 d 
x dx 

o 



d V2 
dx x 
2 

~^ 








1 8.J.4 9 1 
X* dx dx ~x^ 



(18) 



where both operators are manifestly self-adjoint. 

We now make the following crucial observation. It can 
be verified directly that the operator D can be factorized 
as D = —RR , where 



R 



# 

OX 



v/2 








' x 2 dx x 



(19) 



This factorization immediately allows the dynamo equa- 
tions to be transformed into self-adjoint form. Let us in- 
troduce the vector W such that H — EW . As can be di- 
rectly checked, with this definition the vector H automat- 
ically satisfies the solenoidality condition, i.e., its com- 
ponents in the basis I|14ll6fl can be represented as H = 



' d t W 2 ' 




d t W 3 _ 





4c£x 2 

ox 



^9_ r V2 „2_9_S__a_„2 

~ x dx^ x 3 x dx x 4 - dx x 



w 2 
w 3 



,(24) 



and the relation H = EW reads: 



K = --L-^W 3 ). (25, 

Equations (|24|l and (|25|l are the main result of this sec- 
tion. 

For completeness, we note that the equations for the 
functions Mj, and K were first derived by Vainshtein & 
Kichatinov [1J| in the non self-adjoint form: 



1 d 



c 4 dx 



dM 



d t M = —r — x 4 k— + GM - 4hK, 



dx 



(26) 
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Here we adopt the standard notation k — 2i] + kl(0) — 
k l (x), h = g(0) - g(x), G = n" + 4k' /x, and M = M L . 
These equations also follow from H24|) . We also note that 
the Fourier-space versions of i|2tj| I27[) were derived by 
Kulsrud & Anderson for the limit of large magnetic 
Prandtl number (ratio of fluid viscosity to resistivity), 
and by Berger k, Rosner ^3] f° r the general case. 

4. Discussion and conclusion. In systems with no ki- 
netic helicity (i.e., with C(x) = 0), magnetic dynamo 
action is always possible if the magnetic Reynolds num- 
ber is large enough 0, ^| . Since the "helical" terms in 
(|24|l have a destabilizing effect, systems with kinetic he- 
licity should exhibit dynamo action as well. A rigorous 
analysis of the dynamo mechanism requires knowledge 
of the exact spectrum of (|24Jl which is not known for a 
general velocity correlator, k^(x). However, the typical 
behavior of the solution can be understood as follows. 
We introduce the mean-field growth rate, Ao = ffo/ K o- 
From the asymptotic behavior of system l|24|l as x — ► oo, 
one can show that its eigenmodes with A > Ao are local- 
ized, and the closer the growth rate to Ao, the larger the 
correlation length. On the other hand, the eigenmodes 
corresponding to A < Ao have infinite correlation length. 

The formal analogy between Eqs. I|24|) and imaginary- 
time quantum mechanic suggests that the eigenfunctions 
with A > Ao correspond to "particles" trapped by the po- 
tential provided by velocity fluctuations, while the eigen- 
functions with A < Ao correspond to "travelling parti- 
cles". In the non-helical case, where only trapped par- 
ticles have positive eigenvalues, the spacing between the 
eigenvalues decreases with increasing magnetic Reynolds 
number, see e.g., [J^. It is reasonable to expect that the 
same result holds for the helical case. 

We now explain the extent to which the mean-field 
equation (0) describes the dynamo mechanism. Remark- 
ably, in the Kazantsev model, equation (J2J can be derived 
exactly, see, e.g., which allows one to find its pre- 

cise relation to Eq. Ij24(l . In the derivation, B(x,i) is the 
field averaged over the statistical ensemble of the veloc- 
ity fluctuations 10, and the coefficients in the mean- field 
equation (J2J) are given by a — go, f3 = t] + ko/2 ||. In this 
model, A = 270. Eq. (0) can therefore be used to de- 
rive the evolution equation for the correlator of the mean 
field, H l \x,t) = (B l {x,t)&{Q,t)), where the brackets 
denote averaging over the random initial conditions of 
the magnetic field. 

It can be checked that this evolution equation for- 
mally coincides with the large-scale asymptotic (x — > 00) 
of our system (|2"4"j) : consequently it can only be used 
to obtain the large-scale asymptotics of the solutions of 
Eq. (12411 . More precisely, representing the magnetic field 
as B l (x, t) = &(x, t) + SB i (x, t), whereJB* is the fluctu- 
ating part, we can write H lj (x,t) = (B z (x, t)B J (0, t)} + 
(5B i (x,t)5B j (0,t)). We note that while the system 
describes the exact function H i: >(x,t), the mean-field 
equation @ captures only its slowly growing non- 



fluctuating part , which is described by the large-scale 
asymptotics of H i; > (x, t), since the correlation of the fluc- 
tuations vanishes for infinite scale separation. 

In summary, we have used the Kazantsev model to 
compare the exact spectra of magnetic energy and he- 
licity with the predictions of the a-model ©. We have 
demonstrated that the large-scale asymptotics (2 — * 00) 
of the exact solution is described by the non-localized 
eigenmodes (A < Ao) of the self-adjoint dynamo equa- 
tions Ij24(l . This asymptotics can also be derived from the 
mean- field ct-dynamo equation J5J. However, model (J2J) 
misses the faster growing eigenmodes with A > Ao, which 
are present in (|24|l . The correlation lengths of these 
eigenmodes are generally not small. They fill the range 
of scales from the system scale to the the resistive ones, 
so these modes may not be removed by a small-scale 
smoothing procedure. In numerical simulations or as- 
trophysical applications these modes may dominate the 
slowly growing "mean-field" modes. The correct descrip- 
tion of the dynamo mechanism thus requires the resolu- 
tion of the whole range of scales available to the magnetic 
field. 
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